% clear;clc;close all;
% ds = dir('../shots');
% ds = ds(3:end);
% dat = [];
% th = [];
% fh = [];
% dt = 0.004;
% for i = 1:length(ds)
%     [tdat, tth, ~] = ReadSegy(fullfile(ds(i).folder, ds(i).name));
%     dat = [dat tdat];
%     ttrh = [];
%     for j = 1:length(tth)
%         u.offset = tth(j).offset;
%         u.SourceX = tth(j).SourceX;
%         u.GroupX = tth(j).GroupX;
%         ttrh = [ttrh, u];
%     end
%     th = [th ttrh];
% end
load data.mat;